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We consider the problem of characterizing genuine multiparticle entanglement for permutation- 
ally invariant states using the approach of PPT mixtures. We show that the evaluation of this 
necessary biseparability criterion scales polynomially with the number of particles. In practice, it 
can be evaluated easily up to ten qubits and improves existing criteria significantly. Finally, we 
show that our approach solves the problem of characterizing genuine multiparticle entanglement for 
permutationally invariant three-qubit states. 
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I. INTRODUCTION 

Quantum state analysis of large-scale systems is non- 
trivial since many generic methods only work for, at 
most, small system sizes. Experimentally, the manipu- 
lation and control of several qubits has already become 
standard, and current records comprise, for instance, en- 
tanglement between 14 qubits in ion traps [T], ten-qubit 
entanglement using hyper-entangled photons [2] or the 
generation of eight entangled photons [3l 0]. Even for 
such medium scale systems, the available analysis tools 
can be rather cumbersome like, for instance, the task of 
quantum state tomography, i.e., the process to determine 
the underlying quantum state by suitable measurements 
and hence gaining full information. The standard Pauli 
tomography scheme [5] scales exponentially such that a 
feasible, generic tomography protocol for 14 qubits is out 
of scope. 

However, quite often one intends to work with special 
classes of states only. This offers the possibility that one 
can tailor or optimize the analysis tool for those more 
restricted sets. Such more efficient tomography proto- 
cols have been recently designed for generic states of low 
rank [6], particularly important low rank states like ma- 
trix product states [7] or multi-scale entanglement renor- 
malization ansatz states [5], or for states which possess 
some further symmetry like permutation invariance [9]. 

Though full information on a quantum state is appeal- 
ing, it is usually dispensable since one is often more in- 
terested in a few key characteristics or properties of the 
states, which are also used to compare different experi- 
ments. From the plethora of interesting characteristics, 
the main objective of multipartite systems lies on genuine 
multipartite entanglement [TUl |TTJ . This is the strongest 
phenomenon of quantum mechanical correlations within 
such systems that cannot be explained via sufficient con- 
trol on systems of less particle size, known as biseparable 
states. Despite its importance, characterization and de- 
tection of this kind of resource is still hard and only re- 
cently some methods have been proposed [T2TfT7] . A very 
promising detection method constitutes the concept of 
PPT mixtures [T51 US] , the generalization of the positive 



partial transpose (PPT) criterion [2U] to the multipartite 
setting. 

In this paper we tailor, similarly to the tomography 
protocols, the detection of genuine multipartite entan- 
glement via PPT mixtures to permutationally invariant 
states. We show that the question whether a given per- 
mutationally invariant state possesses a PPT mixture 
requires resources which scale only polynomially in the 
number of qubits. Thus, in combination with the to- 
mography protocol [S] (and its variants [HH2H]) and its 
efficient state reconstruction algorithm [21], we develop 
an additional tool to analyze the data after such a quan- 
tum state tomography process. At this point, we would 
like to stress that the derived detection method does not 
rely on the fact that the underlying state indeed pos- 
sesses this symmetry: If the permutationally invariant 
part of a quantum state is entangled, then the complete 
state must be entangled, too [§] • As a further result we 
prove that the criterion of PPT mixtures is necessary 
and sufficient to decide whether a given permutationally 
invariant three-qubit states is genuine multipartite entan- 
gled or not. Thus we obtain another interesting class of 
states, similar to graph-diagonal states of three and four 
qubits [25 , where this approach completely solves the 
question of genuine multiparticle entanglement. As ex- 
amples, we study states like Greenberger-Hornc-Zeilinger 
(GHZ) and Dicke states and obtain strongly improved de- 
tection statements for up to ten qubits. We would like 
to add that in the present paper, we focus on the nu- 
merical evaluation of the criterion of PPT mixtures. Of 
course, also an analytic approach via the construction of 
appropriate witnesses is possible. This can lead to cri- 
teria which can be used for arbitrary particle numbers. 
Results on this problem will be reported elsewhere [2E] ■ 

The structure of this paper is as follows: Section [II] 
summarizes the background on multipartite entangle- 
ment, the concept of PPT mixtures and permutationally 
invariant states. The theoretical results of this paper 
are given in Section III, in particular the aforementioned 
results about structure and scaling of permutationally 
invariant PPT states and about the sufficiency state- 
ment for the three-qubit permutationally invariant case. 
Section |IV| collects some details of our numerical imple- 
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mentation via semidefinite programming (SDP), which is 
used afterwards to test and to compare our method on 
special family of states in Section [V] Finally, we summa- 
rize in Section |VH 

II. PRELIMINARIES 
A. PPT mixtures 

At first, let us review the concept of PPT mixtures [18] . 
which represents a method to detect genuine multipartite 
entanglement. Similarly to the PPT criterion of the bi- 
partite case [2U] , it is a suitable relaxation on the level of 
quantum states which gets more tractable. 

Let us first explain the method for a system of three 
particles because this already highlights the idea. A tri- 
partite state is separable with respect to the bipartition 
A\BC if it can be written as 

Pa\bc = E 9k \ct>1) (41 ® \fl>% a ) & k Bc\ , (1) 
k 

where qk form a probability distribution, i.e., J^k <7fc = 1 
and qk > for all k. The definition is analogous for the 
other possible bipartitions B\AC, C\AB. A biseparable 
state is now defined as a convex combination of states 
which are biseparable with respect to a specific biparti- 
tion, more precisely the state can be written as 

bs sep . sep . sep /o\ 

Pabc = PiP A \bc + PiPb\ac + P*Pc\ab- ( 2 ) 

If a state is not biseparable, it is called genuinely multi- 
partite entangled. 

Naturally, since the bipartite separability problem is 
already hard, the full characterization of biseparable 
states can only be worse. The idea of the PPT mix- 
tures is to define a set of states which includes the set 
of biseparable states, but which is much easier to char- 
acterize than the latter. At the core of this method lies 
the PPT or Peres-Horodecki criterion, which was first 
introduced in Ref. [20]. This criterion, based on the op- 
eration of partial transposition, is a simple and power- 
ful method to detect bipartite entanglement. If a bipar- 
tite system pab is expanded in a chosen tensor prod- 
uct basis as pab — J2ijki PijM (kl\, its partial trans- 
pose with respect to the first subsystem is defined as 
Pab = HtjkiPiiM \kj) {il\. The PPT criterion says that 
if a state pab is separable, its partial transpose is posi- 
tive semidefinite, i.e., it has no negative eigenvalues, or, 
in other words, one says the state is PPT. This implies 
that if has one or more negative eigenvalues (the 
state is then called an NPT state), it must be entangled. 
Throughout this work we will often use the term positive 
meaning positive semidefinite. 

The generalization of this criterion to the multipartite 
case, as introduced in Ref. [18], goes as follows: Simi- 
larly to the definition of a separable state with respect 
to the bipartition A\BC, one can define a state p A P \ BC to 



be PPT with respect to that partition, and similarly for 
the other bipartitions. In analogy to the definition of a 
biseparable state, a PPT mixture of a three party state 
is defined as a convex combination of PPT states with 
respect to a specific bipartition 

pmix ppt . ppt , ppt /o\ 

Pabc = PiPA\bc + PWb\ac + P3Pc\ab- ( 3 ) 

From the PPT criterion, we know that all separable 
states for a fixed bipartition are contained in the set of 
all PPT states for the same bipartition. This then im- 
plies that the set of PPT mixtures contains the set of 
biseparable states. Consequently, if a state is not a PPT 
mixture, it is genuinely multipartite entangled. In Fig.[T] 
we can see a schematic representation of the set of PPT 
mixtures and the set of biseparable states for this three 
particle case. 




FIG. 1: Here we see, for three particle states, a scheme of 
the biseparable states with respect to the three possible bi- 
partitions A\BC, B\AC and C\AB (surrounded by solid blue 
lines). Each of these is contained in the PPT states for each 
respective bipartition (surrounded by dashed green lines). 
The larger solid blue line enveloping the other solid blue lines 
represents the set of biseparable states and similarly the larger 
dashed green line defines the set of PPT mixtures. 

The generalization to the iV-partite case is straightfor- 
ward. A general biseparable state of A*" parties can be 
written as 

P" S = E Pm\mP7w ( 4 ) 

all bipart. 
M\~M 

where the states are separable with respect to partition 
M C {1, . . . , N} and its complement M. The sum runs 
over all possible bipartitions M\M of the N particles. 
Analogously, a PPT mixture of N parties is defined by 

ppmix _ ^ p M — p PPt_ = ^ P m \m, (5) 
M\M M\M 

with operators Pm\m = PMp^m that are positive and 

I A'l | Ai 

PPT with respect to M, 
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The motivation of using the concept of PPT mixture 
is that the problem of determining whether an arbitrary 
state is a PPT mixture or not can be formulated in terms 
of a semidcfinite program (SDP) [57], which makes it an 
easy to implement criterion to detect genuine multipar- 
tite entanglement. This SDP reads as follows 

min — s (6) 
s.t. p = 2J Pm\m> 

M\M 

P M{M > si, P m] m Tm > st, for all M\M, 

where the notation P M ijj — s -"- means that P M m—sl is a 
positive semidchnite matrix. If the result of the optimiza- 
tion s op t is non-negative, then the state p is a PPT mix- 
ture. Otherwise, it is genuinely multipartite entangled. 
Note that in Ref. [18] , the semidcfinite program is written 
in the so-called dual form, which can be interpreted as a 
search for appropriate entanglement witnesses which are 
non-negative on all PPT mixtures. However, both forms 
are equivalent and thus detect the same states. 

In a standard computer, though, this SDP can only be 
applied to generic states of up to five or six qubits [35] . 
In fact, it can be seen that the difficulty of this program 
scales exponentially with the number of qubits. Note 
that the number of different incquivalent bipartitions of 
a system of N particles is given by 2 N ~ 1 - 1. For each 
of these bipartitions the Hermitian matrix P M \jj has 4^ 
free parameters, which are the variables of the semidefi- 
nite program (in fact, one of the matrices is fixed because 
of the equality constraint, but all the others are free). 

We will see in the following that if we restrict ourselves 
to permutationally invariant states, both the number of 
partitions and the number of SDP variables scale only 
polynomially. Furthermore, by an efficient decomposi- 
tion of the matrices which have full or some permutation 
invariance it is possible to check the positivity conditions 
in terms of smaller blocks (see Sec. II B and IV I. This fi- 
nally allows us to construct a SDP able to detect genuine 
multipartite entanglement for larger systems. 



B. Permutationally invariant states 

Many experiments that aim at creating genuine multi- 
partite entanglement, are designed in such a way that the 
generated state is invariant under particle interchange. 
Famous examples are the GHZ or the Dicke states. Math- 
ematically, for any density matrix p a permutationally 
invariant (PI) density matrix can be constructed via 



P(i...n) = [p]pi=w ^OOp^O 



(7) 



where V(ir) is a representation of the permutation tt G 
Sjv acting on the Hilbert space of N qubits. The brackets 
(X...N) should denote invariance under permutations 



between any of the N qubits. We will sometimes employ 
the notation [ . ] pI to explicitly refer to this operation 
in order to shorten the expressions. It can be seen that 
Eq. Q implies that 



[p] PI = V(*)[p] PI V(n)1, for all^eS 



N- 



(8) 



A natural basis to write permutationally invariant states 
is given by the coupled spin basis, for which such states 
attain a particular simple block diagonal form, e.g., [241 
I2§] . In this basis, the Hilbert space of N qubits is de- 
composed as 



H = (C 



2\®N 



N/2 



(9) 



3=3n 



with jmin € {0, 1/2} depending on whether N is even or 
odd. Here, Tlj arc the spin Hilbert spaces of dimension 
2j + 1 and JCj are called the multiplicative spaces, whose 
dimension is given by 



N 



dirmX,)- , N/2 _. 



N 

N/2 



(10) 



for j < N/2 and dim(K. N / 2 ) — 1. The advantage of this 
decomposition is that any permutation V(ir) will only act 
non-trivially onto the multiplicative spaces, that is, they 
can be written as 



N/2 

V{u)= 1®V,-(tt), 



(11) 



where Vj(tt) is an irreducible representation of Sn acting 
on JCj [3171 13T] . This will become important shortly. Fi- 
nally, we will denote the basis states by \j,m,a.j), where 
\j, m) <E Hj and \aj) 6 K-j. These states \j, m, otj) are 
eigenstates of J 2 and J z , where J is the total angular 
momentum operator, while J z is the projection of J in 
the z direction, with 

J 2 \j, to, otj) = h 2 j(j + 1) \j, m, exj) , (12) 
J z \j,m,aj) = hm\j,m,aj) . (13) 

Any permutationally invariant state P(\...n) can m this 
formalism be written as, e.g., [541 129j . 



N/2 



N/2 



P(i... N )= p 3 -pi® dwjh = ® (1 4 ) 

j=j m in J j=j m in 

with states pj of "Hj and a probability distribution pj. 
Complete knowledge of all probabilities Pj and all states 
Pj gives a complete characterization of P(i...n)- Further- 
more, we defined 



B, 



PjPj 
dim(/Cj) ' 



(15) 



which are the blocks that appear in the diagonal of the 
PI state pn...N)- Each of these blocks Bj appears in 
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P(1...N) 



FIG. 2: Due to the form given by Eq. (141, a PI state has a 



block diagonal structure as shown above, in a suitably ordered 
basis. Each block Bj appears dim(/Cj) times in the diagonal. 



the diagonal exactly dhsxQCj) times as shown in Fig. 
From this structure, it is straightforward to compute trie 
number of parameters needed to define a PI state given 
by 



N/2 

E 



(2J + 1) 2 



0=0n 



N + 3 
N 



0(N 3 ), (16) 



which is much smaller than the 4^ — 1 parameters needed 
to characterize a general state. Apart from the better 
scaling, this block structure is also very important for the 
formulation of the SDP. It will be shown in Sec. II VI that 
all constraints of the SDP can be translated into appro- 
priate constraints of the blqcks. For instance, pn...N) > 
is equivalent to Bj > for all j. This is the main reason 
of the polynomial scaling in the end. 



Before we continue we will show how to prove Eq. ( 14 ) 
and, at the same time, a way to carry out the operation 
[ . ] PI without computing the matrices V(w). In general, 
this would be very hard since, for a N particle system, 
V(ir) has dimension A N and there are N\ different per- 
mutations that need to be considered. Although we are 
always dealing with quantum states in the derivation, 
this block structure appears for any permutationally in- 
variant operator. Starting with a general density matrix 



jj'mm' 

(XjOLjl 



(17) 



we have 



[p]pi= E P?/mm' Hi- m ' a ^ 0*'. ft' Wpi 
jj mm 



V 3 {n)\a 3 ){a r \V},{n) 



i i mm' w ' ^ I 



liCjSjj' tr(|oy) (otj>\) 




dim(/Cj) 
j,m) (j,m'\ <E> l/c r 



(18) 



In the first step, we used Eq. (11), while in the second 



we made use of Schur's lemma. Our proof of Eq. (14) 
is similar to what is shown in Ref. [23] except that here 
we obtain explicitly the entries of the blocks Bj from the 
entries of p. In order to see clearly the meaning of the 
result obtained, we define 



B" 3 — > p" 3 "^ , \ j m) (j,m'\ 



such that, this way, the result simply reads 

bV 



dim(/Cj) 



(19) 



(20) 



This means that to calculate the blocks of pn...N) one 
has to take the average over the multiplicative spaces of 
the blocks of p associated with the angular momentum j. 



III. PPT MIXTURES OF PI STATES 
A. Characterization of PI PPT mixtures 

Here, two of the main analytical results of this work are 
presented, in the form of two observations. In the first 
one, we derive a simplified equation which characterizes 
a PI PPT mixture, the PPT mixture of a PI state. We 
show that without losing generality we can restrict the 
sum over bipartitions only to the ones with different num- 
ber of particles on one side. Furthermore, we can impose 
symmetries on the unnormalized PPT states that need 
to be considered. The second observation proves that 
the number of parameters necessary to characterize a PI 
PPT mixture is of 0{N 7 ). 

Observation 1. While Eq. |5]) characterizes a general 
PPT mixture, the equation that characterizes a PI PPT 
mixture can without loss of generality be written as 



N'/2 



pmtx 



P{1...N) 



fe=l ireSjv 



E V(Tr)Q {1 ... mk+1 ... N) V(rf, (21) 



where Q(i...k)\(k+i...N) * s an unnormalized PPT state for 
the partition l...fe|fc+l...JV which is additionally invariant 
under permutations among the first k or the last N — k 
qubits. We have N' £ {N,N — 1} if N is even or odd. 



Proof. Let pn,,,N) be a permutationally invariant state 
which is a PPT mixture. Then, combining Eq. Q and 
Eq.|7| we can write 



pmix 

P(l...N) — 



E Jn E V(n)P MlJI V(n)l (22) 



M\M V(n)eS N 
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Now, let \M\ denote the number of elements inside the 
partition M. Then, for any bipartition M\M, with \M\ = 
k, there is always a permutation tm G 5jv which maps 
M\M to l...k\k + 1...N. We define 



Qi.. 



fc|fc+l...JV 



E V(t m )P miW V(t m )\ (23) 



M\M: 
\M\=k 



which is a positive operator whose partial transpose of 
the qubits 1 ... k is also positive. This is true since 



V(tm)P m] mV(t m ) 



holds for each term in the decomposition and Pm\m being 



PPT for this bipartition. We can now simplify Eq. (22) 
via 



tfL?N) = E ni E V ^ t m t m)P M \m v ^ t m 1t m) 

M\M ' w£S N 

= EjE V(n')V(T M )P MlM V(T M yV(ny 

M | M ' tt'eSn 



N'/2 



E w E ^(^)Oi...*ifc+i...^(^)t 



JV! 

fe=l Tr'GS 



(25) 



where in the first step we defined the permutation n' as 
"kt^I . Furthermore, the Qi..,k\k+i...N can without loss of 
generality be chosen to be invariant under any permuta- 
tion 7Tfe G Sk of the first k or any ttj. € S% of the last 
k = N — k qubits without altering the property that it 
is PPT. This follows since we can actively use the sym- 
metrization as 



E ^MQi...fcifc+i...ivno t 

'n 

fc| E V ( 7T7T k l7T k)Ql...k\k + l...NV(TTTT' k 1 TT k ) 

= h E 1/ ( 7r ')^(^)Qi...fc|fc+i...JV^(7r fc )ty(7r 



1 

~ jfe! 



= E ^ / ( 7r ) < 3(i.../c)|fc+i...Ar^(7r) t , 
where we used the notation 



(26) 



Q(l...k)\k+l...N—y E ^ / ( 7r *:) ( 3l.--fe|fc+l...JV^( 7r fc) t , 



PfcSSfc 



(27) 



and defined the permutation it' as mrk- Note that since 
the partial transpose acts only on the first k particles this 
state is still PPT. Analogously, we can also symmetrize 
Q(i...fc)|fc+i...JV f° r permutations of the last N—k particles 
to obtain Eq. (21). This concludes the proof. □ 



Observation 2. The number of parameters needed to 
characterize a permutationally invariant PPT mixture is 
of 0(N 7 ), where N is the number of qubits. 



Proof. This scaling is due to the fact that Eq. (21) ex- 
hibits two simplifications when compared to Eq. (j5j): The 
first one is that the number of bipartitions that need 
to be considered is only N'/2 and the second is that 
Q(i...k)\(k+i...N) is permutationally invariant within each 
side of its respective bipartition. 



More specifically, if cr^ k ^ and a 



'(k+l...N) 



are operator 

basis elements for permutationally invariant operators of 
k and N—k respectively, any operator Q(i...k)\(k+i...N) 
can be written as 



!(l...k)\(k+l...N) 



E^^a-fe) 

a/3 



a 



(fe+l...JV)- 



(28) 



As mentioned in Sec. |IIB| any permutationally invari- 
ant operator on k qubits can be parametrized by C(fc 3 ) 
parameters. Hence, the operator given by Eq. ( |28[ ) has 
about O [(N — k) 3 k 3 ] parameters, which, at most, can 
be 0(N 6 ) since k can be roughly N/2. This, together 
with the fact that one has to consider about N/2 bi- 
partitions, leads to an overall number of parameters to 
describe a PI PPT mixture of 0(N 7 ). This finishes this 
observation. □ 



B. Necessity and sufficiency for PI three-qubit 

states 

Next we show that the method of PPT mixtures is not 
just only necessary, but also sufficient for biseparability 
of a permutationally invariant three-qubit system. Note 
that this result does not extend to systems of more par- 
ticles where explicit counterexamples are known |32| 133] . 

Observation 3. A permutationally invariant three-qubit 
state is biseparable if and only it is a PPT mixture. 

Proof. Any biseparable state is also a PPT mixture as 



explained in the Sec. II A Thus we are left to show that 
a PPT mixture of a three-qubit permutationally invariant 
state is indeed biseparable. 

For that we can without loss of generality assume the 
special form as given by Observation 1. Since we only 
have non-trivial bipartitions of 1 vs. 2 particles, we ob- 
tain the following form pfjj^?) = [Qa|(bc)]pi or i more 
explicitly, 



p 
Pi. 



(ABC) - g(pA|(BC) + VabPA\(BC)Va 

t 



AB 



VacPa\(bc)Vac)' 



(29) 



where Vab, Vac refer to appropriate permutations. Here 
Pa\(bc) stands for a PPT state with respect to par- 
tition A\BC, which additionally remains invariant un- 
der exchange of system B and C. The structure of 
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permutationally invariant states implies that the two 
qubits BC couple to a spin-1 system, given by the 
symmetric subspace Sym(BC) spanned by |11) , = 
(|01) + 10))/v / 2,|00), and the spin-0 antisymmetric part 
\ip ) = (1 01) - |10))/V2- Thus this state (as well as its 
partial transposition with respect to A) can be decom- 
posed into two parts as 

PA\(BC) = qV A \ Sym (BC) + (1 ~ l)^A ® |VO {lp~ I ■ (30) 

with <J A\Sym(BC) being PPT. Since the symmetric sub- 
space is three dimensional, this state is effectively a qubit- 
qutrit system for which PPT is equivalent to separabil- 
ity [33]. Hence the state of Eq. ( [30| is separable and 
consequently the PI PPT mixture biseparable. □ 

This Observation complements the results of Ref . '2h\ , 
where an analogue result was shown for states with a dif- 
ferent symmetry, namely graph-diagonal states of three 
and four qubits. 



IV. DETAILS OF THE SDP 

Via the simplified form of a PPT mixture of a generic 
PI state the corresponding SDP can now be formulated 
as 



(31) 



s.t. 



P(l...N) 



N'/2 

£ 

k=l 



Q 



a...fc)i(fc+i...jv) 

Ti . . . k 



..k)\(k+l...N)\ pj) 
> Si, 



(l...fc)|(fe+l..JV) 



> si, 



for all k. 



Apart from the polynomial scaling of the number of pa- 
rameters involved in the SDP (c.f. Observation 2), it is 
also crucial to guarantee that the effort needed to ver- 
ify the constraints scales also polynomially. For this, it 
is very important that the parameters are organized in 
blocks for if they were spread throughout the matrix in 
an unstructured way, checking its positivity would still be 
exponentially hard. Thus, in order not to deal with ma- 
trices of size 4 W we want to write the SDP constraints as 
constraints on smaller blocks that constitute the matri- 
ces involved in the program. It can be seen from Eq. ( 28 ) 
that the matrices Qn...k)\(k+i—N)i although not fully per- 
mutationally invariant, must also have a block structure, 
in a suitably chosen basis 



We know from Sec. 



the operator basis elements 
the block structure, 



and 



'(k+l...N) 



1TB] that 
have 



fc/2 



7 (i...fe) 



c« h ®t 



3=3k min 
k/2 

i + i...iv)= D i®^ 

3=3Jimin 



(32) 
(33) 



From Eq. ( 28 ) it then follows 



Q(l...k)\{k+1...N) 

= E^0(^<M^^) 

a,P 3k, 3k 



(34) 



3k, 3k 



where we defined „•_ 

3k, 3k 



C al3 C jk ( 



)D P , . The struc- 

3k 

ture of Eq. (|34| is similar to Eq. (Il4|) which shows that, 
in a suitably ordered basis, the operator Qii...k)\(k+i...N) 
is also block diagonal. Each block Bj k ■ has dimension 
(2jfc + 1) (2jj. + 1) and appears exactly dim(fCj k ) dim(/C JS ) 
times in the main diagonal. In this way, the SDP needs 
to store only the different blocks and keep track of how 
many times each block appears. 

The task now is to define the SDP in terms of the 
these blocks B^ k ^_ which constitute Q(i...k)\(k+i...N) and 
the corresponding blocks Bj of the given PI state pti...m. 
This task is direct for the matrix inequality constraints 
which translate to corresponding matrix inequalities 



r>k 

3k ,3k 



(35) 



for all blocks. Note that the operation of partial trans- 
position is easy to implement, again due to the tensor 
product structure of Eq. ( [28] ) , because 

Qji::.k)\(k+i...N) = J2 Ca 0^t-k)f ® °ffc+i...iv) ( 36 ) 

a p 

where T denotes the usual transposition. 

The implementation of the equality constraint, though, 
requires more care. We know that [Q(i...fe)|(fe+i...jv)] PI IS 
block diagonal in the coupled spin basis of all N par- 
ticles \j,m,oij), as the original PI state, but the ba- 
sis in which each Qn..k)\(k+i...N) i s block diagonal is 
a different one. In fact, it is diagonal in the basis 
\jk,m k ,a 3k ]j^m- k ,a h ) = \j k ,mk,a jk ) ® \j' k ,mt,a jf ). 
In general, the basis transformation between two spins 
jk 7 jk to a combined total spin j is given by the Clebsch- 
Gordan coefficients 



li, m) 



-3k<m k <jk 
~3k< m k<3k 



(jk , m k ; jk , m h |i, m) \ j k , m k ; j h , m h ) 



.(37) 

This transformation holds for each spin, irrespective of 
the multiplicative spaces. However, via the multiplicative 
spaces one keeps track of how many spins jk of system 
1 . . . k, and how many spins j k of system k + 1 . . . N cou- 
ple with each other. To compute the resulting blocks 



B 



3 0f Efe 

follows: 



[Qd. 



.k)\(k+l...N)_ 



] the procedure is hence as 



From the blocks of Q(\ 



fe)|(fe+i...iv), typically 
given in the basis \j k , m k \ j k , mj) , one first com- 
putes their contribution to each total spin j, using 
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\j,m), via the transformation of Eq. (37). This re- 
sult is denoted as B-'^_. Note that one only gets 
a non-trivial matrix if the two individual spins can 
form at all a total spin j, i.e., \jk~3k\ — 3 — 3k+3k- 



• Afterwards one performs the average over all pos- 
sibilities, more precisely, 



B 3 = dim (fcj J dim (fcjfc) g fc,j 



jkjk 



dim(/Cj 



(38) 



Here dim(/Cj fc ) dim(/Cj- ) is the number of spins 
3k, jk (which couple to a spin j) in the original op- 
erator. 

This way, after summing over k, the RHS of the equality 
constraint of Eq. (31) is computed. Due to its symme- 
try, Q(i...fe)i(fe+i...jv) has a polynomial number of param- 
eters, so the basis transformation requires the compu- 
tation of only a polynomial number of Clebsch-Gordan 
coefficients. This discussion leads to the final observation 
of this work: 

Observation 4. The SDP to detect genuine multipartite 
entanglement of PI states of N qubits via the concept 
of PPT mixtures can be formulated in terms of 0(N 3 ) 
matrices whose size is at most 0(N 2 ). 

Since the number of different blocks of a general PI op- 
erator on k parties is 0(k), each operator Qn...k)\(k+i...N) 
has about 0(N 2 ) different blocks. Since we need roughly 
N/2 of these operators for a PI PPT mixture we have of 
0(N 3 ) operators Bj k j_ in total. For each of these oper- 
ators we need to check two matrix inequalities. Further- 
more, the biggest of the blocks is the one with k — N' /2 
and j k = N'/i,^ = (2N - N')/4, hence the maximal 
dimension is 0(N 2 ). 



EXAMPLES 



White noise tolerance: N-qubit Dicke states with k excitations |D N k > 
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Huber et al. 
PPT mix 



FIG. 3: Comparison between the PPT mixture criterion 
(filled symbols) and the criterion of Huber et al. [13: (empty 
symbols) of the white noise tolerance for Dicke states \Djsr,k)- 
We show the white noise tolerances for N up to 10 qubits 
and k up to N/2. The PPT mixture criterion is more ro- 
bust to white noise and in some cases, the difference of the 
white noise tolerance between both criteria is very significant 
reaching values larger than 40%. 



Example 1. — Dicke states have first been studied in 
the context of light emission from a cloud of atoms [35] 
and have been prepared in many experiments |36l I37j . 
The symmetric -/V-qubit Dicke state with k excitations 
is defined as the superposition of all basis states with k 
excitations, 



\D 



N.k 



® |0) 



®N-k 



(39) 



In this section, we present some examples for the ap- 
plication of the SDP, illustrating the strength of the PPT 
mixtures. Our first example is the calculation of the 
white noise tolerance for Dicke states which we compare 
to the method of Ref. [13]. In our second example, we 
consider a mixtures of a GHZ state, a W state and white 
noise and compare the detection range with the ones of 
Refs. [12 E] , which were the best known results so far 
for this class of states. In both cases, we achieve sig- 
nificantly improved results. However, since our method 
is based on a numerical approach, it is limited by the 
memory of the computer and we could only run it for 
states of at most 10 qubits with the first prototype of 
the program. Note that we did not further optimize the 
algorithm for these special kind of states. In contrast the 
criteria of Rcfs. [12j [14] are analytic and can therefore 
be applied to arbitrary qubit numbers. Further criteria 
from the approach of PPT mixtures, which similarly rely 
on analytic estimates will be discussed elsewhere [2S] , 



and is therefore a permutationally invariant state. For 
example, the four-qubit Dicke state with two excitations 
is given by |£> 4 , 2 ) = (|0011) + |0101) + |0110) + |1001) + 
1 1010) + |1100>)/V6. 

We computed the white noise tolerance for Dicke states 
of up to 10 qubits and up to N/2 excitations comparing 
it to the criterion of Ref. [13] ■ From Fig. 3 we see that 
the PPT mixture criterion is always more robust to noise 
and the difference is more significant for larger number 
of qubits and excitations. The improvement reaches val- 
ues larger than 40%, which should prove itself useful for 
current Dicke state experiments. 

Example 2. — Other well known states, which are also 
invariant under permutations are the GHZ and W state. 
The GHZ state is defined as 

\GHZ N ) = j= (\0f N + \lf N ), (40) 





PPT mix 




Refs. [12,14] -- 


Biseparable 





0.2 0.4 0.6 0.8 1 

P1 

FIG. 4: Here, we show which three-qubit states (lower trian- 
gular) as defined by Eq. ( |42[ ) are detected by the PPT mix- 
tures criterion as genuinely multipartite entangled (area out- 
side the solid line). We compare it to the criteria presented in 
Refs. |12lll4j . consisting of two inequalities, one optimized for 
the GHZ state and the other for W state (dashed lines). The 
PPT mixtures criterion represents a significant improvement 
and is optimal for three-qubit PI states, so the states inside 
the solid line are biseparable. 




FIG. 5: The same as Fig. [4] but for eight-qubit states. The 
PPT mixture criterion represents a very significant improve- 
ment, but note that here we cannot prove that the PPT mix- 
ture criterion is necessary and sufficient for biseparability. 



VI. CONCLUSIONS 



while the W state is the Dicke state with one excitation, 
\W N ) = ■— : (|10...0) + |01...0) + ... + |0...01)) . (41) 

In this example, we consider the following ./V-qubit states 

1 

P[Pl,P2) =PlPGHZ N +P2PW N + (1 -Pi -P2)^v, (42) 

with PghZn = \GHZ N ){GHZ N \, PWn = \W N )(W N \. 
These states are a convex combination of a GHZ state, a 
W state and white noise. 

Such states can be represented by a point in a two- 
dimensional plane whose coordinates are given by pi and 
P2, as shown in Fig. [4] and Fig. [5} Of course, only certain 
pairs of pi,p2 correspond to valid quantum states, hence, 
only the lower triangle shown in these figures describe ac- 
tual quantum states of this class. In these figures, which 
corresponds to three and eight qubits respectively, we 
furthermore show the set of states which is detected by 
the PPT mixture criterion in comparison to the method 
developed in Refs. [TUQ3]. I n both cases the set of states 
detected by the PPT mixtures is much larger than the 
one verified by the aforementioned criteria and this dif- 
ference grows with the number of qubits. In Fig. [4} in 
fact, the PPT mixtures criterion is optimal, cf. Sec. |III B] 
so the states that have a PPT mixture are biseparable. 



In this paper we tailored the detection of genuine mul- 
tipartite entanglement via PPT mixtures for permuta- 
tionally invariant states. In contrast to a generic ./V-qubit 
state, where the question of characterizing PPT mixtures 
scales exponentially with the number of particles, our op- 
timization for this special class of states only requires a 
polynomial scaling of the order 0(N 7 ). This was possi- 
ble by deriving a more restricted but still general form 
of a PPT mixture using the additional symmetry of the 
state. Via this method, we were able to analyze more 
rigorously the entanglement for system sizes where the 
original numerical PPT mixture method would fail. 

In addition, we have shown that the criterion of PPT 
mixtures completely solves the question of genuine multi- 
partite entanglement for permutationally invariant three- 
qubit states. This furthermore supports the conjecture, 
motivated in Ref. [25] . that PPT mixtures are necessary 
and sufficient for biseparability of three qubits. We leave 
this open for further discussion. 

On more general ground, we believe that the develop- 
ment or optimization of existing analysis tools to larger- 
scale systems is a mandatory step for a well-grounded 
investigation of the properties of systems with many par- 
ticles and its experimental implementation. This should 
help to close the gap between methods for small system 
playgrounds and the really interesting system sizes that 
could deserve the term quantum computer at some time. 
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